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Qv^ ■ Gamma-ray bursts may occur within pulsar wind bubbles (PWBs) under a 

o\ : 
■ 

\o ■ 
o : 

. it comes to death and produces a fireball. We here study the early afterglow 

o ■ 

emission from GRBs expanding into such a PWB environment. Different from 
Qh the usual cold GRB external medium, the PWBs consist of a hot electron- 

positron (e + e~) medium with typical 'thermal' Lorentz factor of the order 
c/3 . of j w , the Lorentz factor of the pulsar particle wind. After GRB blast waves 

r* ■ shock these hot e + e pairs, they will emit synchrotron radiation peaking at 

GeV bands. It is shown that GcV photons suffer negligible absorption by 
H ' the soft photons radiation field in PWBs. Thus, strong GeV emissions in the 

early afterglow phases are expected, providing a plausible explanation for the 
long-duration GeV emission from GRB940217 detected by EGRET. Future 
GLAST may have the potential to test this GRB- PWB interaction model. 
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number of scenarios, such as the supranova-like models in which the progenitor 
pulsar drives a powerful wind shocking against the ambient medium before 



1 INTRODUCTION 

Although the fireball shock scenario of GRBs and their afterglows has been well demon- 
strated by multiwavelength observations of afterglows (e.g. Wijers, Rees & Meszaros 1997; 
Piran 1999; van Paradijs, Kouveliotou & Wijers 2000) , the nature of the central engine is not 
yet well known. There is observational evidence for the association of long duration bursts 
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with star-forming regions (Bloom et al. 2002), or possibly with supernovae (SNe; Galama et 
al. 1998). One class of candidates of the central engine involves that the supernova leaves 
behind a rotationally-supported supermassive or massive neutron star, which then shrinks 
by shedding angular momentum via a pulsar-type wind and finally collapses to a black hole 
surrounded by a disk (Supranova model; Vietri & stellar 1998) or transits to a strange star Q 
(Wang et al. 2000), producing the GRB fireball. The pulsar- wind bubble (PWB) forms when 
the relativistic wind from the pulsar shocks against its ambient matter and creates a pulsar 
nebula like the famous Crab nebula. Though different physical processes are involved for 
the death of the neutron star in these two scenarios, the PWB and the SN ejecta behaviors 
as well as the rotation energy of the progenitor pulsar are nearly identical; hence the two 
scenarios have rather similar implications to the evolution of GRB afterglows. This class 
of models [j] in which SN explosion preceded the GRB event has recently gained support 
from the detections of strong Fe emission features in the X-ray afterglow spectra of some 
GRBs (e.g. Piro. et al. 1999;), particularly that of GRB991216 (Piro et al. 2000), as well 
as an absorption feature in the prompt emission of GRB990705 (Amati et al. 2000). These 
detections implies that a large amount (<; O.1M ) of pure iron is located in the vicinity 
(r ^ 10 16 cm) of the GRB source (e.g. Vietri et al. 2001). Such a surprising large iron mass 
is most naturally produced in a SN explosion and the inferred distance of the SN ejecta 
implies the SN event precedes the GRB by months to years or shorter (see, e.g., Bottcher, 
Fryer & Dermer 2002). 

Effects of wind nebulae resulted from the central 'pulsar' on the GRB afterglows have 
recently been examined by Inoue et al. (2002) and Konigl & Granot (2002; hereafter KG) 
from different points of view. Assuming the SN ejecta may be fragmented by the plerion 
radiation, Inoue et al. (2002) study the inverse Compton effect of this radiation field on 
GRB external shock decelerated by the outlying baryonic material. Different from Inoue et 
al. (2002), in the model of KG, the afterglow shock wave propagates in a bubble composed 
mainly of shocked electrons and positrons that emanate from the central pulsar. KG argue 
that the pulsar wind nebula environment, bounded by the SNR shell, can account for the 

* Strange quark matter is conjectured to be more stable than hadronic matter (Witten 1984). The existence of strange matter 
is allowable within uncertainties inherent in a strong-interaction calculation (Farhi & Jaffe 1984), Strange stars, composed of 
this kind of quark matter, may exist and could be born from a massive neutron star as it spins down (Wang et al. 2000) 
t In the literature, Fe emission features are usually interpreted in the supranova model (Vietri & Stellar 1998). We note however 
that Fe line production is irrelevant to whether the progenitor pulsar collapses to a black hole or converts to a strange star. 



© 2002 RAS, MNRAS 000, 000-000 



Early GeV afterglows from gamma-ray bursts in pulsar wind bubbles 3 

high electron and magnetic energy fractions (e e and e#, respectively) inferred in a number 
of afterglow sources, as the PWBs are expected to have a significant e + e~ component and 
to be highly magnetized. Following the previous treatment of the PWB structure that was 
developed for the Crab and other nebulae (e.g. Rees & Gunn 1974; Kennel & Coroniti 1984, 
hereafter KC84), and assuming the magnetic pressure remains in approximate equipartition 
with the particle pressure throughout the shocked-wind bubble, they developed the structure 
model of PWB. On the other hand, a one- zone simplified version of the shock-wind model 
of PWB was developed by Chevalier (2000) and is shown to be well consistent with the X- 
ray luminosity of nebulae. Guided by these works, we attempt to give a simplified analytic 
description of the PWB that may exist in the supranove-like models (section 2), and study 
the early afterglow emission from blast wave propagating in the hot e + e~ bubble (section 3) 
. The characteristic synchrotron frequencies of early afterglows are much higher than that of 
the afterglows in a typical interstellar medium (ISM). In section 4, we give a brief summary. 

2 PULSAR WIND BUBBLES 

Before collapsing, the supermassive or massive neutron star loses its energy through a pulsar- 
type wind, whose luminosity can be estimated from the magnetic dipole formula L w = 
B 2 R%Q 4 /6c 3 = 3 x lO^ergs- 1 Bf 2 R 6 eq l3 Qj and the spin-down time, in which half of the 
neutron star rotation energy AE rot ~ 10 53 erg need to be lost, is 

tsd = AE rot /L w = 10 yr A£ , rotj53 S^ 2 i2~ jl3 fiJ , (1) 

where B = 10 12 -Bi 2 Gauss, R eq = 13i? eg) i 3 Km and f2 = 10 4 Sl 4 s~ 1 are, respectively, the dipole 
magnetic field, the typical equatorial radius and the angular velocity of the neutron star. Due 
to the unknown variations of these parameter values, we below consider t sd as a free param- 
eter and parameterize the wind luminosity as L w = AE rot /t s d = 3 x 10 44 ergs~ 1 A_E roti53 t^ 1 , 
where t sd = lOyrt^i. 

We assume that the PWB is elongated along the rotation axis of the pulsar with an outer 
radius Rb equal to that of the SN ejecta in this direction and pulsar wind is shocked at radius 
R s (Rees & Gunn 1974; KC84; KG). The pressure force of the expanding PWB is expected to 
accelerate the SN ejecta into a thin shell (Reynolds & Chevalier 1984) moving with a velocity 
v e j = 0.05— 0.1c as a significant amount of AE rot is deposited into the bubble. Approximately, 
the outer radius of the nebula at time t sd equals to Rb = Vbt S( i = 10 18 cm 1^0.1^,1 , where 
Vb = v e j = O.lci^o.i is the expansion velocity of the outer edge of the PWB. Fe line emission 
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requires that the pre-ejected SN ejecta should be no farther than r ~ 10 16 cm. To reconcile 
this short Fe line-emitting distance with the radial distance of the PWB, KG suggest that 
the PWB become elongated in the polar direction because of anisotropic mass outflow from 
the GRB progenitor star. We note that the anisotropic pulsar wind power could also give 
rise to the elongated configuration. In fact, some pulsar wind nebulae have been found to 
show such anisotropic structure( e.g. Gaensler et al. 2001). 

The pulsar rotation energy is assumed to go into a highly relativistic wind of the Lorentz 
factor 7 W . Spectral modelling of the Crab nebula has given a current value of ~ 3 x 10 6 
for 7k,. However, people do not know the initial value of r y w for pulsars in their early years. 
A much lower initial value 7^ < 10 4 for Crab pulsar is inferred by Atoyan (1999) when 
interpreting the radio spectrum of the Crab nebula. In light of these factors, we regard that 
the value of 7^ is quite uncertain and may span a range of two orders. The wind luminosity 
consists of electromagnetic and particle contributions with their flux ratio denoted by the 
magnetization parameter a. KC84 note that a must be smaller than 0.1 in order to have 
a high efficiency of the wind power into synchrotron radiation. We assume that the pulsar 
wind flux is dominated by e + e~ pairs, so L w ~ 47m 1 7 2 i? 2 m e c 3 , where n\ is the proper pair 
density in the wind just before the shock front. 

As the young SNR expands at a velocity Vb <C c, the wind energy excluding the part that 
goes into radiative and adiabatic loss will accumulate within a volume confined by SNR. The 
characteristic radius R s of the shock front is determined by the balance between the ram 
pressure of the wind L w /4ttcR^ and the total magnetic and particle pressure P in the bubble. 
Since the sound speed in the shocked fluid c/ is much larger than that of the SN ejecta, no 
pressure gradient in the nebular bubble can be maintained for dynamically important time, 
and the bubble is virtually isobaric (Reynolds & Chevalier 1984). At the shock front R s , the 
wind e + e~ particles may acquire a power-law energy distribution of the form ^(7) oc 7~ p 
for 7 > 7 m , where n — J N( / -f)d , -f is the number density of the shocked relativistic particles, 
7 is the Lorentz factor of the particles and 7 m is the minimum value of 7. Usually p > 2 
is assumed to keep the total energy of the shocked particles finite and p = 2.2 is inferred 
for Crab nebula from its X-ray spectrum (e.g. Pravdo & Serlemitsos 1981). p ~ 2.2 seems 
also consistent with the simulation result of the first-order Fermi particle acceleration (Kirk 
et al. 2000). From the shock jump conditions, the internal pressure immediately behind the 
shock front is pb = nim e c 2 7^, and the energy density is e& = 3nim e c 2 7^. As in Chevalier 
(2000), we assume that the energy density in the downstream region is divided between a 
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fraction e e into particles and a fraction eg = 1 — e e into the magnetic field, similar to the 
definition in the blast- wave model of GRB afterglows (Sari et al. 1998). The plausibility for 
incorporating the equipartition assumption have been accounted for in Chevalier (2000) and 
KG. 

With these assumptions, the magnetic field in the bubble is 

B b = (^) 1/2 = 0.17 Gauss AE 1 ^^ (2) 

where a = R s /Rb and e# = 0. 5^0. 5. The minimum Lorentz factor of the shocked e + e~ 
is 7m = (p - 2)/(p - l)e e 7™ = 8 x 10 3 e ei0 . 5 7™,5, where e e = 0.5e ei0 . 5 and 7™ = 10 5 7 Wj5 . 
The synchrotron lifetime of the shocked e + e~ with the average Lorentz factor % = e e ^( m is 
estimated to be 

247TT77<g C g 2 2. 1 13 a 1 12 / \ 

Tsyn = a T l {B b ) 2 = X a (Ee > - 5eB > - 5 ^ w <^ sd > 1 ^^rot,bZ V bfi.li (3) 

where Ot is the Thompson cross section. Thus, for t sd ^ lOyr and 7^ ranging from 10 4 to 
10 6 , the synchrotron lifetime of the e + e~ in the bubble is much shorter than the age of the 
bubble and the PWB is in the strongly cooling regime for parameters of interest. In this case, 
the radial width of the bubble A b is much smaller than the radius of the bubble R b , which 
means a ~ 1, and we can express the bubble volume as V = AttR 2 A b . Now let's estimate 
the radial width A b at the time (t s d) when GRB-PWB interaction begins. Equating the ram 
pressure with the internal pressure of the bubble, we have 
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where e p is the internal energy density of the particles. The particle energy density can be 
regarded as the sum of the contributions from both the e + e~ that have suffered strongly 
cooling and those that have not suffered any cooling yet, i.e. 

e P =( e e L w ^dt+ e e L w dt)/4%R 2 b A b , (5) 

W0 7 e Jt sd -T syn J 

where 7^ is the Lorentz factor of cooled e + e~ in the bubble at time t s d- 7^ = 67rm e c/ (aT{B b ) 2 (t sd — 
t)). Finally we get 

mZ £ \ ln T^ + I' 3 X ^"^KS.s^as^i^^io.i. ( 6 ) 

MT\£> ) tBlw T syn / <J 

consistent with our previous assumption A b <C R b . 

The innermost part of the bubble is composed of newly injected pairs from the pulsar 
that have suffered negligible cooling and the minimum Lorentz factor is 7m - a nonmag- 
netic ultrarelativistic shock, the downstream fluid just behind the shock front moves with 
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a velocity ~ |c relative to the shock front. The bulk of the shocked e + e~ may move at 
a speed between and the bubble expanding speed v b . As a rough estimate, the width of 

these 'uncooled' e+e" region is A x ~ vr syn ~ 10 16 cm7- 1 5 e~o. 5 esl).5^d,i A£; r _ t,53^,o.i' where 
we have taken v being 0.2c. 

3 GRB-PWB INTERACTION AND EARLY GEV AFTERGLOWS 
3.1 Early GeV afterglows and GRB940217 

At the time t s d after the SN explosion, the pulsar collapses to a black hole or converts to 
a strange star and the GRB goes off, sending a fireball and a relativistic blast wave into 
the PWB. Different from the usual cold proton-electron preshock medium in the standard 
model of afterglows, the PWB is hot, composed of randomly moving, relativistic e + e~ pairs, 
and highly magnetized. The initially freely-expanding fireball shell will be decelerated by 
the swept-up hot e + e _ pairs. We can calculate the extent to which the shell is decelerated by 
PWB gas. The total enthalpy, representing the effective inertia that decelerates the shell, is 
U & 4irR 2 b A b (e b + Pb ) ~ 1.2 x 10 51 erg^^ ge e g 5 e s ^Q 5 tg rf l A£/ r0 ( 53. If we temporarily neglect the 
radiative loss of the GRB shell (i.e. the shell energy is conserved), its Lorentz factor will be 
reduced from the initial value T to T f = {E^/Ufl 2 ~ lOiEfa) 1 ^^.^ J^AE^ 
when the shell collides with the outer supernova remnant, where Eq = 10 53 erg_Eo 53 is the 
initial isotropic kinetic energy of the GRB shell. The deceleration of the GRB shell starts 
when the energy of the swept-up medium equals to Eq/Tq, lagging the GRB prompt emission 
by a time Ati ag = R b /2Tlc = 180sec ^1^0.1X0,300) where r = 300Po,3oo is the initial 
Lorentz factor. As the forward shock expanding into the hot PWB, the average Lorentz 
factor per GRB-shocked e + e~ particle is % = e e 7„,r = 1.5 x 10 7 e ei0 .57 w ,5r3oo ; if the post 
GRB-shocked medium has the same equipartition factors e e , e# as those in the pulsar wind 
shocks. According to the shock jump conditions, the minimum Lorentz factor of the GRB- 
shocked e + e~ is 

7™ = j—^ e e 7»r = 2.5 x lO^o^^rsoo (7) 

for p = 2.2 and the magnetic field in the post GRB-shocked medium is 

B s = YB b = SIGauss r 30 o^L^f A^^o.i- (8) 

The shock- accelerated electrons may have a maximum Lorentz factor 7^ determined by the 
balance between the shock accelerate time (which could be the gyroperiod) and the radiative 
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loss time. Taking the shock-acceleration time as the gyroperiod and assuming synchrotron 
cooling , we get 

37Tg e 1/a _ 7 -1/2 -1/4 3/4 ar-1/4 1/2 ^ 

IM—y^gs) — iy X iU 1 300 e B,0.5 I .sd,l^ A - C/ roi,53 t 'b,0.1; l J J 

where g e is the electron charge. The observed characteristic synchrotron emission photon 
energy for e + e~ with the Lorentz factor 7^ is 

K. = T ^ 2 ^~c = 1 - 1GeVr 30oT^5ee 1 o.5eKL^ 2 A<5 5 3^o.i ) (10) 

and the maximum synchrotron emission photon energy is 
In 2 

hv s M = -^^T = 45GeVr 30 o, (11) 
m e coT 

with the spectrum being F v oc v~ v l 2 between them ^|. 

When the GRB blast wave swept-up matter cools, an initial GeV burst occurs at the 
deceleration radius. The deceleration of the GRB shell starts when the energy of the swept- 
up medium equals to Eq/Tq, and the corresponding radial distance that the shell travelled 
in the PWB is Xdec = (Eg/TlAnRKeb + Pb)) = 3 x 10 13 cmro| o*sd,i- Please note that, as the 
shell slows down in the PWB, the radial distance x that the shell travelled is much shorter 
compared to its radius. In this case, the angular spreading time scale [j] of the shell emission 
At ang = Rb/2T 2 c (Katz 1994a; Fenimore et al. 1996) dominates over the radial time scale 
^krad — x/2T 2 c, Thus an interesting situation arise: though the time that the shell spends 
on moving through the PWB is very short, the shock generated energy radiates over a much 
longer observed time which is determined by the shell angular spreading time. Therefore the 
shell angular spreading time is relevant to the flux and light curve in the following calculation. 
So, the effective deceleration time should therefore be td ec = Rtf^^c = 180secfb ) o.it s d,iPo~3oo- 
As the shocked PWB medium is in the fast cooling regime, the power emitted is simply 

t Note that for a range of parameter values, the assumed 7J may become larger than jm> then one may ask what will happens. 
To answer this, one must know the detailed microphysics of how the shock is formed. If there is a process which thermalises 
particles faster than the gyroperiod, then this process could lead to higher energy 'thermalised' particles rather than the 
'nonthermal' maximum produced by the standard Fermi acceleration picture (J. G. Kirk 2001, private communication). In this 
case, the Fermi picture would not be relevant. If there is no such a process, then the Lorentz factor of the GRB-shocked e+e — 
pairs implied by the jump conditions cannot be reached. The energy must be radiated away within the shock structure itself 
(by synchrotron radiation) and the jump conditions must be modified to account for the energy 'leak' (J. G. Kirk 2001, private 
communication). Thus, strong high-energy emissions are also expected. 

3 Because of relativistic beaming, an observer sees up to solid angle of T _1 from the line of sight. Two photons emitted at the 
same time and the same radius from the shell, one on the line of sight and the other at an angle of T^ 1 away, travel different 
distance to the observer. The difference leads to a delay in the arrival time by At ang = Rf,/2F 2 c . 
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that given to the shocked e + e~, that is e e times the power generated by the shock, i.e. 
J F v dv = e e ^§-4^T where d^ is the luminosity distance of the GRB source. Integrating the 
fast cooling spectrum with v s m v s c over frequency (Sari et al. 1998) , we have 

dE s 1 



/ 



2 + -?-ri-(2kr* 



p-r~ w m J J 6e dt Airdr (12) 

where E s is the isotropic kinetic energy of the decelerated GRB shell. Strictly, to calculate 
the radiated power at the deceleration radius, we need to calculate the total number of 
swept-up relativistic electrons and their radiative power as in Dermer et al. 1999 . But as 
a rough estimate, we assume that the bulk of the shell kinetic energy is dissipated into the 
shock energy at the deceleration radius (Rees & Meszaros 1992; Dermer et al. 1999) and get 
the fluence at v s m ~ lGeV band (forp = 2.2): t dec vF u {lGe\[) = 8 x 10" 6 ergcm~ 2 e ei0 .5-E! 3 G^ 2 28 
. Clearly, this intensity of GeV emission is well above the detectability of EGRET detector. 
This synchrotron GeV emission will fade with time and to a certain time it will be too weak 
to be detected, leading to an extended GeV afterglow emission. 

Below we will study the fading behavior of GeV afterglow. According to Eq. (12), the 

observed afterglow flux at ~ lGeV is 

v dF s 
vF u (lGeV) = v s m FL(-r p/2+1 oc W 2 " 1 ^- (13) 

The duration of emissions from the shell scales with r as t ~ Rb/2T 2 c oc T~ 2 . As the 
shocked gas is partially radiative, the shock energy E s decreases as E s oc T e in the limit 
that r 3> T 3> 1 (Bottcher & Dermer 2000), where e is the fraction of the shock-generated 
thermal energy that is radiated. For a fast-cooling shock, e = e e . Considering v s m oc T 4 , from 
Eq.(13) we find the GeV spectral power flux decays with time as 

^(lGeV) ~ 4.4 x 10-Wm- 2 s- 1 (^)-^ 1 )-^/ 2 ee ,o. 5 ^3<28Q!i^o 1 iro,30o- (14) 

For comparison with EGRET observation and future GLAST observation, we now cal- 
culate the energy fluence the detector can collect during a certain time duration t (see Fig. 4 
of Zhang & Meszaros 2001), i.e. the product of vF v Gev flux and the observing time. 

ti/F„(lGeV) ~ 8 x 10- 6 ergcm- 2 e e ,o.5^3<28(T^) _( ^ 2K£e/2 - (15) 

''dec 

The fluence threshold for EGRET is roughly ~ 2.1 x 10 _6 ergcm -2 for short integration time 
regime (Zhang & Meszaros 2001). According to Eq. (15), the decreasing of the fluence from 
the initial value at time td ec to the sensitivity threshold of EGRET implies that the observed 
GeV emission lasts 
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+ ~ f ?c;nn c ^ t 2/(a4+0 - 5ee ' a5) l'P s \2/(0.4+0.5e e ,o. S ) J-4/(0-4+0.5e e , . 5 ) . p-2 /-, fi N 

for p = 2.2. We think that this mechanism provides a plausible explanation for the famous 
long-duration (~ 5400 sec) high-energy photons, with energies ranging from 36 Mev to 18 
GeV, from GRB940217 detected by EGRET (Hurley et al. 1994). The future GLAST de- 
tector is much more sensitive than EGRET and its fluence threshold is roughly ~ 1.2 x 
10 _9 t 1,/2 ergcm -2 for a long integration time regime, so a more extended detected GeV emis- 
sion time is expected with t dur ~ 1.8 x 10W^ 4+a5ee '°- 5) ^^ 

^6,o.i^sd,iro~3oo f° r V — 2-2 . The power-law fading behavior of GeV emissions in our scenario 
distinguishes itself from some other models, such as afterglow inverse-Compton model ( 
Meszaros & Rees 1994; Dermer 2000; Zhang & Meszaros 2001) and hadron process models 
(e.g. Katz 1994b) 

3.2 Synchrotron self-Compton and nebula-induced external Compton 
components at GeV band 

We need to consider the synchrotron self-Compton (SSC) and nebula-induced external 
Compton (EC) components at GeV band in order to compare it with synchrotron com- 
ponent obtained in the above subsection. The GRB-shocked e + e~ gas will inverse-Compton 
scatter synchrotron photons emitted by themselves to higher energies. Because the shock is 
fast-cooling, the ratio of the SSC to synchrotron luminosity is approximately given by (Sari 
& Esin 2001) 



Lssc 



B > (17) 

^ 1 L syn if e e > e B 



For fast-cooling shocks, L syn ~ ^nd\v s m Fl m and L S sc ~ knd\ u^ l sc F^ c , respectively. Note 
that ^> rn e c 2 , the SSC energy emission peaks at u KN = Y^ s M m e c 2 

~ 3TeVr3QQe^ 1 // 5t^ 4 L A£^ 7 ^^ 5 4 3 t'^' 2 1 , above which the Klein-Nishina effect suppresses Compton 
scattering. Therefore, at ~ lGeV band, the intensity of the SSC fluence component is 
generally below that of the synchrotron one so long as e e is not much larger than e#, though 
at TeV band, it dominates over the latter. 

Another important process that might contribute to the high-energy component is the 
Comptonization of external nebular soft photons by relativistic e + e~ in the expanding GRB 
shell. The luminosity of the EC emission can be obtained analogously to the SSC case, by 
using ~ T^Ydm^m^ura where r e is the optical depth of the shocked material. The 
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number density of e + e~ in the bubble is n b e = 37„,ni = 3L W / ' 47i r y w R 2 m e c 3 = 0.037"^^ A£' rot i;^' _ 1 , 
so r e ~ l/3cr T n b e Ax ~ 10~ w ^~\ l tJ^ 1 AE rot v^Q A Axi 6 , where Ax is the radial distance the shell 
travelled in the bubble. The luminosity of the PWB v h m F b m can be obtained from Eq. (9) 
in Chevalier (2000): u b m F b m ~ lO-^ergcm^s^ei^.^^QV^ for p = 2.2. We finally get 
the EC luminosity at time td ec , L EC ~ 1.1 x 10 _12 ergcm~ 2 s _1 , for typical parameters used 
. Noting that T^O^JJ 3> m e c 2 , the EC energy emission peak also at v KN = Y^ s M m e (? ~ 
3Te"VT3Qoe^^ \t^ x AE~^^ z v]j^ x . It is clear that EC contribution to ~ lGeV band is also 
unimportant compared to the synchrotron component. 

3.3 Attenuation of high-energy photons by soft photons in PWB 

As the PWB has plenty of soft photons, we must consider whether the high-energy photons 
suffer strong attenuation due to pair-production reaction while propagating through the 
PWB. The optical depth of a high-energy photon with energy E 1 due to 7 — 7 absorption 
writes r 77 = f^™ ax (r(v>)n(i>)Abdi> , where hv = 2{m e c 2 ) 2 /E 1 and u max is the frequency of 
the maximum energy photon in the bubble, n(v) ~ L{y) /hvAnRfc is the number density of 
photons with energy hu (L(is) is the luminosity per unit frequency of the PWB), and a(u) 
is the pair production cross section of E 1 with soft photon hu: a{y) ~ (T t(^) _1 f° r v > v . 
The formula of the PWB luminosity L(u) has been derived by Chevalier (2000) (see his Eq. 
(9) ), given by 



/ , p-2y- 1 


6g e 2 


(p 


\p-l) 


.(2nm e c) 2 c_ 





^-l c g-2)/4 7£ _2 i? -(p-2)/2 L J + 2)/4 1/ - P /2 j (lg) 

for v > v b m > v\ and L{y) oc h>~ 1 / 2 for v h m > v > v b . Setting r 77 = 1, after a simple algebraic 
computation, we get hu Q ~ O.lSeV^eiS^T^^^A^^^oT for p = 2.2 . Therefore 
the high-energy cutoffs due to pair production with soft photons in PWB lie at 

TP ^ 7-^-0.18^0.86 „,0.73 + 0.14 a iti-0.045^,0.09 /iq\ 

E cut ~ 51eVe e 5 e B)0 . 5 7 w , 5 ^,i AA r 0t ,53 v b,o.i ( 19 ) 

for p = 2.2. So we can safely conclude that the synchrotron GeV photons can reach us without 
significant attenuation due to pair production through interaction with soft photons in the 
PWB. 



4 SUMMARY 

GRBs could arise within PWBs in supranova-like scenarios (Vietri & Stella 1998; Wang et 
al. 2000) of GRBs, in which SN explosion precedes GRB by a time from several months to 
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tens of years and the pulsar emits a relativistic wind, making a wind bubble while shocking 
against the ambient medium. After the spin-down time, the central pulsar may collapse 
to a black hole or convert to a strange star, forming the GRB fireball. In this paper, we 
studied the early afterglow high-energy emissions from GRB blast wave expanding into 
such a pulsar wind bubble. We find that, owing to that the PWB is composed of hot 
e + e~ pairs with high randomly-moving Lorentz factors, the majority of the GRB-shocked 
PWB pairs have such high thermal Lorentz factors that their characteristic synchrotron 
frequencies lie at GeV bands. Our calculation shows that the attenuation of GeV photons by 
soft photons in PWB due to pair-production reaction is negligible. Strong power-law decaying 
GeV emissions during the early afterglow phase is naturally expected in this scenario and 
it may provide a plausible explanation for the long-time (~ 5400s) GeV photons from 
GRB940217. The predicted power-law decaying behavior of the flux is distinguished from the 
afterglow electron-IC mechanism (Zhang & Meszaros 2001) as the latter predicts a delayed 
hump in the light curves. 
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